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SUMMARY 

Interpolation  techniques  provide  an  alternative  to  the  phase-locked 
frequency  aultlplier  for  the  calculation  of  the  tine  doaain  average  of  a 
vibration  signal.  Higher  order  interpolation  techniques  produce  flatter 
passbands  and  lover  sidelobes  in  the  stopband  but  require  longer 
calculation  tines.  Aliasing  errors  are  introduced  into  the  result  by 
replication  of  the  sidelobes  during  interpolation.  In  general,  the  errors 
are  attenuated  by  tine  doeain  averaging,  but  under  soee  conditions  nay  be 
passed  without  attenuation. 
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1.  INTRODUCTION 


Tiae  doaain  averaging  is  a  signal  processing  technique  which  per aits 
the  extraction  of  periodic  vaveforas  froa  noisy  signals.  It  is 
particularly  useful  for  the  analysis  of  the  vibration  of  aechsnical  systeas 
such  as  gearboxes,  as  it  enables  the  vibration  of  a  single  gear  to  be 
separated  froa  the  vibration  of  the  coaplete  systea.  As  the  vibration  of 
the  gear  is  shown  in  the  tine  doaain  over  one  coaplete  revolution, 
differences  in  the  perforaance  of  the  teeth  around  the  gear  becoae 
apparent,  enabling  faults  such  as  fatigue  cracks  to  be  identified. 
Recently,  tine  doaain  averaging  has  been  applied  with  considerable  success 
to  the  detection  of  fatigue  cracks  in  the  input  spiral  bevel  pinion  in  the 
aain  rotor  gearbox  of  the  Vessex  helicopter  [1). 

Tiae  doaain  averaging  differs  froa  the  conventional  technique  of 
spectral  analysis  in  that  it  requires  a  rotational  reference  signal  to 
enable  the  angular  position  of  the  desired  gear  in  the  systea  to  be 
deterained.  In  its  siaplest  fora,  this  reference  signal  aay  consist  of  a 
pulse  train  which  is  synchronized  with  the  rotation  of  the  required  gear. 
This  pulse  train  may  be  used  to  control  the  saapling  of  the  vibration 
signal.  Hovever,  the  acquisition  of  a  signal  directly  froa  the  required 
gear  or  its  shaft  is  often  not  practicable  as  these  coaponents  aay  not  be 
readily  accessible.  In  such  cases,  it  is  necessary  to  take  a  reference 
signal  froa  soae  other  location  in  the  systea  which  is  readily  accessible, 
and  to  convert  that  signal  to  the  required  sampling  control  signal  using  a 
phase-locked  frequency  aultiplier.  For  the  analysis  of  the  vibration  of 
the  Wessex  input  pinion,  a  simple  phase-locked  frequency  aultiplier  was 
developed  to  convert  the  alternator  output  vavefora,  at  a  noainal  400  Hs, 
into  a  pulse  train  which  provides  256  samples  per  revolution  of  the  pinion 
121- 


But  the  phase-locked  frequency  aultiplier  is  not  without  its  own 
problems.  In  the  phase-locked  frequency  aultiplier  developed  for  the 
Wessex  pinion,  an  edge- triggered  digital  phase  detector  was  used.  It  has 
been  found  on  soae  occasions  that  this  is  susceptible  to  false  triggering 
by  electrical  noise  on  the  alternator  signal  or  dropouts  on  the  tape 
recording,  causing  the  loop  to  lose  lock  and  so  corrupt  the  tine  doaain 
average.  Furthermore,  in  order  to  optimise  the  performance  of  the 
phase-locked  frequency  aultiplier  for  any  given  application,  it  is 
necessary  to  select  the  loop  filter  characteristics  to  suit  the  expected 
input  frequency  and  the  aultiplication  ratio  of  the  loop.  Changing  the 
multiplication  ratio  to  enable  the  tiae  doaain  averaging  of  different  gears 
in  the  systea  requires  the  changing  of  the  loop  filter  and,  for  a 
multi-stage  aultiplier,  this  adds  greatly  to  the  complexity  of  the  design. 
Furthermore,  the  phase-locked  loop  requires  a  finite  tine  to  acquire  and 
lock  onto  the  input  signal.  Under  soae  circumstances  this  time  can  be  very 
large,  necessitating  long  recordings  of  the  data  of  vhlch  only  a  small  part 
can  actually  be  used  when  the  loop  has  finally  locked.  Even  when  locked, 
the  loop  must  inevitably  lag  behind  any  changes  in  the  phase  or  frequency 
of  the  input  signal,  and  these  errors  oust  limit  the  accuracy  of  the 
average  which  is  obtained. 


As  an  alternative  to  the  phase-locked  frequency  multiplier,  a 
technique  for  time  doaain  averaging  based  on  off-line  interpolation  by  a 
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digital  computer  can  ba  aaployad.  If  both  tha  vibration  signal  and  tha 
rotational  rafaranca  signal  ara  saaplad  slaultanaoosly  at  a  flxad  clock 
frequency,  than  for  aach  saapla  of  tha  vibration  signal  tha  eorraspondlng 
rotational  rafaranca  position  will  also  ba  known.  Proa  tha  numbers  of 
taath  on  tha  gaars  in  tha  systaa,  tha  rafaranca  signal  can  ba  ralatad  to 
tha  rotation  of  tha  required  gaar,  and  tha  vibration  signal  can  than  ba 
intarpolatad  to  obtain  tha  vibration  of  tha  systan  at  incranonts  of  tha 
rotation  of  tha  raqulrad  gaar.  From  thasa  valuas  tha  tina  donaln  avaraga 
nay  ba  calculatad. 

This  tachniqua  has  savaral  advantagas  ovar  tha  phasa-lockad  fraquancy 
multiplier.  Firstly,  tha  sansitlvity  to  nolsa  on  tha  rafaranca  signal  is 
raducad,  as  it  is  possibla  to  procass  tha  saaplad  reference  signal  to 
remove  sons  of  tha  Interference.  As  all  of  the  processing  is  performed 
off-line,  aore  tine  is  available  for  such  operations  than  if  on-line 
averaging  is  attempted.  Secondly,  only  one  complete  cycle  of  tha  rafaranca 
signal  is  required  before  tha  technique  can  start  producing  Intarpolatad 
data,  unlike  tha  phasa-lockad  fraquancy  aultiplier  which  nay  require  aany 
cycles  before  lock  can  ba  achieved.  Thus  aany  aore  averages  can  be 
obtained  froa  the  saaa  quantity  of  data.  but  the  biggest  advantage 
proaisas  to  ba  tha  ease  of  changing  the  multiplication  ratio.  Vith  tha 
phase-locked  frequency  aultiplier  this  required  the  changing  of  loop 
filters,  but  with  the  interpolation  technique  it  can  be  achieved  in  a 
coaputer  program  simply  by  changing  tne  ratio  relating  the  rotation  of  the 
required  gear  to  the  reference  signal.  In  principle,  this  makes  possible 
the  calculation  of  the  tine  doaain  average  of  any  gear  in  a  mechanical 
system. 

It  is  believed  that  a  linear  interpolation  technique  nay  already  have 
been  used  for  tha  tine  doaain  averaging  of  gear  vibration.  Linear 
interpolation  is  vary  simple  to  implement,  but  the  results  it  produces  are 
only  approxiaata.  It  is  vary  important  that  interpolation  techniques,  as 
they  may  ba  applied  to  tine  doaain  averaging,  should  be  understood  froa  a 
digital  signal  processing  viewpoint,  in  order  that  the  frequency  response 
of  the  techniques  nay  be  estimated  and  that  parameters  may  be  selected  to 
optimize  the  performance  of  the  techniques.  This  report  examines  the 
theory  of  interpolation  as  applied  to  tine  domain  averaging  and  describes 
the  sample  and  hold,  linear  and  cubic  interpolation  techniques,  comparing 
the  frequency  response  of  each.  It  is  revealed  that  interpolation 
techniques  introduce  aliasing  errors  which  produce  images  of  spectral 
components  at  other  frequencies.  It  is  shown  that  while  linear 
Interpolation  can  ba  satisfactory,  cubic  Interpolation  provides  greater 
accuracy  as  it  has  a  batter  frequency  response  and  reduced  aliasing  errors, 
although  it  does  require  a  longer  calculation  time.  The  techniques  are 
illustrated  by  examples  relevant  to  the  analysis  of  the  vibration  of  the 
input  spiral  bevel  pinion  in  the  main  rotor  gearbox  of  the  Vessex 
helicopter. 


7.  FREQUENCY  RESPONSES  OF  INTERPOLATION  TECHNIQUES 


Saapla  and  hold  interpolation,  also  known  as  zero  order  hold,  is  tha 
simplest  interpolation  technique.  The  incoming  signal  x(t)  is  saaplad  and 
tha  value  obtained  is  held  until  tha  next  saapla  has  bean  taken,  after 
which  that  new  value  is  held,  and  so  on.  To  evaluate  the  fraquancy 
response  of  tha  saapla  and  hold  interpolation,  tha  technique  is  expressed 
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in  t«ru  of  simple  digital  aignal  procaaaing  operations.  Piratly,  tha 
incoming  aignal  ia  aaaplad  at  a  frequency  f,  by  multiplying  tha  aignal  by  a 
train  a(t)  of  unit  lapulaaa  with  pariod  Ta-l/f(,  aa  illuatratad  in  Plgura 
la.  By  tha  convolution  theoraa  [3],  aultlplication  in  tha  tiaa  domain  ia 
equivalent  to  convolution  in  the  fraquancy  doaain.  Thua  in  tha  fraquancy 
domain,  tha  Pourlar  tr ana form  X(f)  of  tha  aignal  ia  convolvad  with  tha 
Pouriar  t ran* form  S(f)  of  tha  train  of  impulaaa,  which  is  alao  a  train  of 
lapulaaa  but  at  aultiplaa  of  tha  sampling  fraquancy  fs.  This  producas 
raplicatas  of  tha  original  basaband  at  aultiplaa  of  tha  sampling  fraquancy, 
as  shown  in  Pigura  lb.  Tha  rasult  is  than  convolvad  vith  a  ractangla 
function  a(t)  of  width  T*  and  amplitude  1,  as  illuatratad  in  Plgura  2a.  By 
tha  convolution  thaoraa  (3),  convolution  of  signals  in  tha  time  domain  is 
equivalent  in  the  fraquancy  domain  to  tha  multiplication  of  tha  Pourlar 
transforms  of  thoaa  signals,  producing  tha  rasult  shown  in  Plgura  2b.  Vhan 
intarpolation  is  parforaed  and  the  Fourier  transform  A(f)  of  tha  ractangla 
function  la  multiplied  by  the  Fourier  transform  of  the  sampled  signal,  the 
replicates  must  be  included  in  the  latter.  Parts  of  tha  raplicatas  My  be 
passed  by  the  sldelobes  in  the  Fourier  transform  of  the  rectangle  function, 
although  at  reduced  amplitude. 

The  rectangular  function  a(t)  may  be  defined  by: 


a(t) 


1  for  0  <  t  <  Ts 
V)  for  t  •  0,  t  -  T* 
0  elsevhere 


(1) 


Its  Fourier  transform  A(f)  is  given  by: 


-J2nf(T»/2) 

A(f)  -  (sinfifl^fj/itf )  e 

-J2llf(T,/2) 

-  T,  (sin( nTj f )/( ifTj f ) )  e 


(2) 


Not*  the  presence  in  tha  above  equation  of  the  sidelobes  which  decrease 
inversely  vith  the  frequency.  Note  also  the  non- zero  phase  characteristic, 
which  is  due  to  the  holding  of  the  previous  sample  value,  thereby 
introducing  a  tin*  delay  of  T»/2  into  the  result.  A  simpler  presentation 
of  interpolation  has  been  presented  (4),  but  it  is  incorrect  as  it  neglects 
the  replication  in  the  frequency  domain  caused  by  the  sampling  of  the 
signal . 

Linear  interpolation,  also  known  as  first  older  hold,  is  another 
simple  interpolation  technique.  After  sampling  the  signal,  the  required 
Intermediate  values  are  obtained  by  linear  interpolation  betveen  the 
samples  either  side  of  the  required  abscissa.  To  evaluate  the  frequency 
response  of  linear  interpolation,  the  technique  is  expressed  in  terms  of 
simple  digital  signal  processing  operations.  Firstly,  the  incoming  signal 
is  sampled  as  for  the  sample  and  hold  Interpolation.  Then  the  sampled 
signal  is  convolvad  vith  a  triangular  function  b(t)  of  vidth  2TS  and 
amplitude  1,  as  shown  in  Figure  3a.  In  the  frequency  domain,  the  Fourier 
transform  of  the  sampled  signal  is  multiplied  by  the  Fourier  transform  B(f) 
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of  the  triangular  function,  as  illustrated  in  Figure  3b.  Parts  of  the 
replicates  of  the  baseband  aay  be  passed  by  the  sidelobes  of  the  latter,  as 
for  the  saaple  and  hold  interpolation,  but  as  the  sidelobes  of  the 
triangular  function  are  smaller  than  those  of  the  rectangular  function,  the 
amount  passed  should  be  smaller. 

The  triangular  function  b(t)  may  be  defined  byt 


b(t)  -  1  -  t/Ts  for  |t|  <  Ts 

-  0  elsewhere 

Its  Fourier  transform  B(f)  is  given  by: 

B(f)  -  sin2(itTsf)/(nf),Ts 

-  Tj  (sin<itTsf)/(nTjf))2 


(3) 


(4) 


Note  the  presence  in  the  above  equation  of  the  sidelobes  which  decrease 
inversely  with  the  square  of  the  frequency,  more  rapidly  than  for  the 
rectangular  function.  Unlike  the  saaple  and  hold  Interpolation,  the  phase 
characteristic  is  zero,  as  no  time  delay  is  introduced. 

For  better  performance,  a  more  sophisticated  interpolation  technique 
such  as  cubic  interpolation,  also  known  as  third  order  hold,  is  a  logical 
choice.  Vhereas  linear  interpolation  utilized  a  single  sample  either  side 
of  the  required  abscissa,  cubic  Interpolation  uses  two  samples  either  side. 
A  cubic  function  may  be  fitted  simply  to  four  samples  by  the  Lagrangian 
interpolation  formula  of  order  four  [5J.  To  evaluate  the  frequency 
response,  the  interpolation  is  expressed  in  terms  of  simple  digital  signal 
processing  operations.  Firstly,  the  incoming  signal  is  sampled  as  for  the 
sample  and  hold  and  linear  interpolations.  Then  the  sampled  signal  is 
convolved  with  the  impulse  response  c(t)  of  the  cubic  function,  having  a 
width  4TS  and  amplitude  1,  as  shown  in  Figure  4a.  In  the  frequency  domain, 
the  Fourier  transform  of  the  sampled  signal  is  multiplied  by  the  Fourier 
transform  C(f)  of  the  cubic  function,  as  illustrated  in  Figure  4b.  Parts 
of  the  replicates  of  the  baseband  may  be  passed  by  the  sidelobes  of  the 
cubic  function,  but  as  these  are  very  much  smaller  than  those  of  either  the 
sample  and  hold  or  linear  interpolation  functions,  the  signal  passed  by  the 
sidelobes  will  be  correspondingly  smaller.  Unlike  the  saaple  and  hold  and 
linear  interpolation  techniques  for  which  the  determination  of  the  impulse 
response  was  easy,  for  cubic  interpolation  it  is  rather  more  difficult,  and 
the  impulse  response  is  different  depending  on  where  the  four  samples  are 
placed  relative  to  the  required  abscissa.  Only  the  case  in  which  two 
samples  are  located  either  side  of  the  required  value  will  be  considered 
here.  Techniques  for  calculating  the  impulse  response  of  Lagrangian 
interpolation  functions  have  already  been  described  (6).  Using  these 
techniques,  the  impulse  response  c(t)  for  Lagrangian  interpolation  of  order 


four  is  found  to  be: 

c(t)  -  H2<t/Ts 

♦  2) 

for 

-2TS 

<  t  <  -T, 

=  Mt/Ts 

*  1) 

for 

-Ts 

<  t  <  0 

tk 
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Ho(t/Ts) 

for 

0  <  t  <  Ts 

H-i (t/Ts  -  1) 

for 

Ts  <  t  <  2T; 

where 

H2(v) 
Hx(v) 
H0(v) 
H-i (v) 


(v  +  l)(v)(v  -  l)/6 
-(v  +  l)(v)(v  -  2)/2 
(v  +  1 ) ( v  -  l)(v  -  2)/2 
-(v)(v  -  l)(v  -  2)/6 


The  Fourier  transform  C(f)  is  found  to  be: 

C(f)  =  Ts  (1  +  (2nTsf)2/6)  (sin(nTsf)/(nTsf))4 


(5) 


(6) 


(7) 


Note  in  the  above  equation  the  presence  of  the  sidelobes  which,  for  large 
frequencies,  decrease  inversely  with  the  square  (not  the  fourth  power)  of 
the  frequency.  The  sidelobes  will  be  smaller  than  for  the  linear 
interpolation  due  to  the  fourth  power  of  the  sine  term.  The  phase 
characteristic  is  zero  as  no  time  delay  is  introduced. 

The  frequency  responses  of  the  sample  and  hold,  linear  and  cubic 
interpolation  techniques  are  compared  in  Figure  5.  Increasing  the  order  of 
the  interpolation  technique  decreases  the  amplitude  of  the  sidelobes,  so 
producing  better  stopband  performance.  The  cubic  function  produces  a 
flatter  passband  than  either  the  sample  and  hold  or  linear  functions,  and 
also  has  the  sharpest  cutoff  rate.  Thus  the  cubic  function  clearly 
provides  better  frequency  response,  although  the  higher  order  must  increase 
the  calculation  time. 


3.  ALIASING  ERRORS  PRODUCED  BY  INTERPOLATION 


In  the  previous  section,  interpolation  techniques  were  introduced  and 
their  effects  in  the  time  and  frequency  domains  were  demonstrated. 
However,  the  results  produced  were  continuous  functions  and  not  discrete 
samples,  even  though  the  original  signal  had  been  sampled.  The 
interpolated  results  obviously  cannot  be  used  in  this  form,  as  a  digital 
computer  can  only  operate  on  discrete  values.  These  could  be  created  by 
sampling  the  interpolated  signals,  although  in  practice  this  would  not 
normally  be  done,  as  the  interpolation  function  need  be  evaluated  only  at 
the  required  abscissa,  thus  reducing  the  time  required  for  computation. 
However,  to  simplify  the  analysis,  the  evaluation  will  be  modelled  as  a 
sampling  process. 

The  interpolated  signal  is  sampled  as  shown  in  Figure  6a  by 
multiplying  it  by  a  train  u(t)  of  unit  impulses  with  a  period  Tj ,  where  in 
general  Tx/Ts.  From  the  convolution  theorem,  multiplication  in  the  time 
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domain  is  equivalent  to  convolution  in  the  frequency  domain.  Thus  the 
Fourier  transform  of  the  interpolated  signal  is  convolved  with  the  Fourier 
transform  U(f)  of  the  train  of  impulses,  which  is  also  a  train  of  impulses 
at  multiples  of  the  interpolation  frequency  f^l/Tj,  as  illustrated  in 
Figure  6b.  A  particularly  important  results  becomes  apparent.  The 
sidelobes  which  are  present  in  the  continuous,  interpolated  signal  are 
replicated  during  the  convolution  and  may  overlap  the  baseband  and 
stopband,  so  producing  errors  in  the  result.  These  errors  vill  be  referred 
to  here  as  aliasing  errors,  because  the  image  of  a  signal  is  made  to  appear 
at  a  different  frequency.  These  aliasing  errors  should  not  be  confused 
with  the  more  common  use  of  the  term  to  describe  errors  which  are  produced 
by  the  sampling  of  a  signal  at  an  insufficient  rate. 

The  magnitude  of  the  aliasing  errors  depends  on  three  factors. 
Firstly,  it  depends  on  the  magnitude  of  the  sidelobes  present  in  the 
frequency  response  of  the  interpolation  function,  and  secondly  the 
separation  of  the  replicates,  which  is  determined  by  f t ,  and  thirdly  by  the 
bandwidth  of  the  original  signal.  The  smaller  the  sidelobes  and  the  signal 
bandwidth  and  the  larger  the  interpolation  frequency  then  the  lower  will  be 
the  interference  between  the  sidelobes  and  the  main  lobe.  Reducing  the 
sidelobes  requires  the  use  of  more  sophisticated  interpolation  techniques. 
Increasing  the  interpolation  frequency  requires  more  frequent  evaluation  of 
the  interpolation  function.  Both  vill  increase  the  time  taken  to  calculate 
the  result.  Reducing  the  bandwidth  limits  the  usefulness  of  the  result.  A 
compromise  solution  is  required  which  combines  acceptable  levels  of 
aliasing  errors  with  acceptable  calculation  time. 

In  order  to  demonstrate  the  effects  of  aliasing  errors,  an  example 
will  be  presented  which  is  relevant  to  the  analysis  of  the  vibration  of  the 
input  spiral  bevel  pinion  in  the  Wessex  helicopter  main  rotor  gearbox.  The 
pinion  rotates  at  a  nominal  frequency  of  43  Hz,  while  the  alternator 
signal,  nominally  400  Hz,  is  used  as  the  rotational  reference  signal.  The 
frequency  of  the  gear  rotation  fR  is  related  to  the  alternator  frequency  fA 
by  the  numbers  of  teeth  on  the  gears  and  the  number  of  poles  in  the 
alternator.  After  elimination  of  common  factors,  this  is  given  by  [2]s 


f*  =■  (19/61)  x  (21/61)  fA 


(8) 


The  first  example  features  test  signals  synthesized  as  follows.  The 
alternator  signal  is  modelled  as  a  sinusoid  of  unit  amplitude  and  frequency 
fA-401.57  Hz.  The  vibration  signal  is  the  sum  of  128  sinusoids  of  unit 
amplitude  at  frequencies  from  1  to  128  times  fR  in  increments  of  fR.  Both 
are  evaluated  at  intervals  corresponding  to  a  sampling  rate  fs=20480  Hz  and 
stored  in  blocks  of  512  sample  pairs  on  disc.  The  value  of  fA  was  chosen 
to  produce  exactly  51  samples  per  simulated  alternator  cycle.  Hence  for 
the  test  signal: 


fs  -  51  fA 


(9) 


Applying  the  sample  and  hold  interpolation  technique  to  the  test  signal  and 
sampling  the  result  at  a  frequency  fi=1024fR  produces  a  signal  having  the 
amplitude  spectrum  shown  in  Figure  7.  Many  aliasing  errors  are  visible, 
with  the  largest  of  these  occurring  at  about  350  orders  due  to  the  first 
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upper  sidelobe  of  the  baseband.  Other  major  features  can  also  be 
identified  as  sidelobes  of  replicates.  Some  even  overlap  the  baseband  and 
interfere  with  signal  in  the  band  of  interest  between  1  and  128  orders. 
Using  linear  interpolation  of  the  same  data,  the  spectrum  in  Figure  8  is 
produced,  with  greatly  reduced  aliasing  errors  in  both  the  passband  and 
stopband.  Using  cubic  interpolation,  the  spectrum  in  Figure  9  is  produced, 
with  further  reductions  in  aliasing  errors. 


4.  EFFECT  OF  TIME  DOMAIN  AVERAGING  ON  ALIASING  ERRORS 


In  the  previous  section,  the  effects  of  sampling  the  interpolated 
signals  were  demonstrated.  A  key  feature  was  the  introduction  of  aliasing 
errors  by  the  replication  of  sidelobes  falling  vithin  the  passband  and 
stopband.  For  sample  and  hold  and  linear  interpolation,  these  errors  were 
significant.  Fortunately,  with  many  signals  time  domain  averaging  can  help 
to  reduce  these  errors.  Averaging  is  performed  using  the  interpolated 
signal  evaluated  at  fi  where  fi  is  an  integer  multiple  of  the  rotation 
frequency  fR  of  the  required  gear.  That  is,  fi=LfR  where  L  is  an  integer. 
The  original  sampling  frequency  fs  is  selected  from  the  sampling  rates 
which  are  available  on  the  data  acquisition  system,  and  as  this  is  in 
general  independent  of  the  rotation  frequency  of  the  mechanical  system,  fs 
and  fi  will  in  general  be  independent.  All  of  the  components  passed  by  the 
time  domain  averaging  process  will  be  at  multiples  of  the  rotation 
frequency  fR,  and  as  f:=LfR,  all  of  the  components  in  the  replicates  of  the 
baseband  will  be  at  multiples  of  fB .  However,  aliased  images  of  these, 
present  in  all  of  the  sidelobes,  will  not  in  general  be  at  multiples  of  fH. 

For  example,  consider  the  first  sidelobes  of  the  baseband  produced  by 
the  sampling  of  the  original  signal  at  fs>  A  component  in  the  baseband  at 
a  frequency  fM-MfR,  where  M  is  an  integer,  will  be  passed  by  time  domain 
averaging.  Its  aliased  image  in  the  first  sidelobe  will  appear  at  fs+ffi, 
which  will  not  in  general  be  an  integer  multiple  of  f „ ,  and  so  after 
averaging  the  image  will  be  attenuated.  Similarly,  for  the  first  lower 
sideband  of  the  first  replicate  of  the  baseband,  the  aliased  image  will  be 
located  at  fi-fs-fii,  which  in  general  will  not  be  an  integer  multiple  of 
fR,  and  so  with  averaging  will  be  attenuated.  The  effects  of  averaging  on 
aliasing  errors  will  be  demonstrated  using  the  test  signal  described  in  the 
previous  section.  The  amplitude  spectrum  of  the  signal  obtained  after  128 
averages  using  sample  and  hold  interpolation  is  shown  in  Figure  lOThere  is 
a  marked  improvement  over  that  shown  in  Figure  7  as  the  passband  is  now 
smooth  and  the  aliasing  errors  in  the  stopband  have  been  attenuated.  For 
linear  interpolation,  the  result  is  seen  in  Figure  11,  and  similarly  for 
cubic  interpolation  in  Figure  123. 

In  general,  time  domain  averaging  tends  to  attenuate  aliasing  errors 
which  are  generated  by  interpolation,  but  there  are  specific  instances  in 
which  aliasing  errors  will  not  be  reduced.  For  example,  consider  an 
unwanted  periodic  noise  at  a  frequency  fN,  such  as  might  be  generated  by 
another  gear  meshing  frequency,  with  fN*NfR  where  N  is  an  integer.  Thus  in 
the  baseband,  fR  will  be  attenuated.  However,  sampling  of  the  signal  at  a 
frequency  fs  produces  replication  of  the  Fourier  transform  about  fs»  so 
that  an  image  of  fN  appears  at  fs±fN.  After  averaging,  the  image  at  fs-fH 
will  remain  if: 
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fs  -  ”  K£r  where  K  is  an  integer  (10) 


If  K*350,  then  from  Equations  8,9  and  10,  fN=125.62fR.  Thus  a  periodic 
noise  signal  occurring  at  125.62  orders  will  produce  an  aliased  image  at 
350  orders  which  will  be  attenuated  by  the  filter  shape  of  the 
interpolation  function,  but  will  not  be  attenuated  further  by  time  domain 
averaging. 

This  is  illustrated  by  comparing  the  spectra  obtained  after  1  and  128 
averages  of  a  test  signal  synthesized  featuring  a  periodic  component  at  a 
frequency  of  125.62  orders.  Vith  sample  and  hold  interpolation,  the 
spectrum  after  a  single  average  in  Figure  13  shows  many  images,  with  the 
biggest  consisting  of  a  peak  near  126  orders  surrounded  by  many  sidebands, 
but  vith  a  clean  peak  at  350  orders.  After  128  averages,  the  peak  at  350 
orders  remains  unchanged  in  amplitude  while  all  other  components  are 
attenuated,  as  shown  in  Figure  14.  For  linear  interpolation,  after  a 
single  average  there  are  fewer  images  as  illustrated  in  Figure  15,  with  the 
largest  peak  occurring  again  near  126  orders  and  a  clean  peak  at  350  orders 
having  a  lover  amplitude  than  before  due  to  the  smaller  sidelobes.  After 
128  averages,  the  peak  at  350  orders  remains  unchanged  in  amplitude  while 
all  other  components  are  attenuated,  as  shown  in  Figure  16.  For  cubic 
interpolation,  after  a  single  average  there  are  even  fever  images  as 
illustrated  in  Figure  17.  After  128  averages,  the  peak  at  350  orders 
remains  unchanged  as  shown  in  Figure  18.  This  demonstrates  the  importance 
of  the  amplitude  of  the  sidelobes  in  suppressing  this  class  of  aliasing 
error. 

Errors  can  also  be  produced  in  the  passband.  For  example,  images  can 
be  produced  in  the  second  lover  sidelobe  of  the  replicate  about  fi  such 
that: 


fx  -  2fs  +  fN  =  Kf„  vhere  K  is  an  integer  (11) 


If  K»120,  then  from  Equations  8,9  and  11  vith  fi=1024fR,  fH=47.33fR.  Thus 
a  periodic  noise  component  at  47.33  orders  vill  produce  an  image  at  120 
orders.  The  image  vill  be  attenuated  by  the  filter  shape  of  the 
interpolation  function  but  vill  not  be  attenuated  further  by  averaging. 

This  is  demonstrated  for  1  and  128  averages  for  a  test  signal 
synthesized  vith  a  periodic  signal  at  47.33  orders.  With  sample  and  hold 
interpolation,  after  a  single  average  many  images  are  produced  as  shown  in 
Figure  19.  The  biggest  peak  occurs  near  47  orders  surrounded  by  many 
sidebands,  but  vith  a  clean  peak  at  120  orders.  After  128  averages,  as 
shown  in  Figure  20,  the  peak  at  120  orders  remains  unchanged  in  amplitude, 
while  the  components  near  47  orders  are  heavily  attenuated.  With  linear 
interpolation,  after  a  single  average  no  images  can  be  seen  on  the  original 
scale,  only  the  major  peak  near  47  orders,  as  illustrated  in  Figure  21. 
After  128  averages,  Figure  22  shows  on  an  enlarged  scale  that  there  is  a 
small  component  present  at  120  orders,  together  vith  the  remains  of  the 
peak  near  47  orders  which  has  not  yet  been  completely  attenuated.  For 
cubic  interpolation,  no  images  can  be  seen  on  the  original  scale  in  Figure 
23,  only  the  major  peak  near  47  orders.  After  128  averages,  Figure  24 
shovs  on  an  enlarged  scale  that  there  is  a  very  small  component  present  at 
120  orders,  together  with  the  remains  of  the  peak  near  47  orders.  This 
demonstrates  that  aliasing  errors  can  be  produced  in  the  passband  which  may 
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not  be  completely  attenuated  by  time  domain  averaging.  The  more 
sophisticated  interpolation  schemes,  having  better  stopband  performance, 
help  to  reduce  the  amplitude  of  these  errors. 

The  examples  considered  so  far  have  assumed  that  the  rotation 
frequency  fR  is  exactly  constant.  This  is  unlikely  to  occur  in  practice, 
particularly  in  helicopter  gearboxes  for  which  speed  variations  of  up  to 
±1  X  might  be  expected  at  a  nominally  constant  torque  setting.  These  small 
speed  or  frequency  variations  will  be  referred  to  as  jitter.  Synchronous 
noise  at  a  frequency  fN  will  vary  directly  with  variations  in  the  rotation 
reference  frequency  fA.  The  interpolation  frequency  fx  will  also  vary 
directly  with  fR  since  fT=LfR,  but  the  sampling  frequency  fs  remains 
constant  as  it  is  fixed  by  the  data  acquisition  equipment.  Thus  the 
equations  developed  previously  to  describe  the  images  may  not  hold.  That 
is,  if  fN=125.62fR,  the  aliased  image  will  not  remain  steady  at  350fR  as  fR 
and  fN  vary.  Similarly,  if  fN=47.33fR,  the  aliased  image  will  not  remain 
fixed  at  120fR.  As  a  result,  the  noise  will  be  attenuated  by  time  domain 
averaging.  The  extent  of  the  attenuation  will  depend  on  the  magnitude  of 
the  jitter.  Although  the  worst  case  amplitude  of  the  aliasing  errors  may 
be  reduced,  other  periodic  components  may  produce  images  for  a  short  time, 
thus  contributing  to  low  level  noise  in  the  average.  It  is  difficult  to 
predict  the  amplitude  of  the  aliasing  errors  without  a  knowledge  of  the 
possible  periodic  components  and  the  frequency  distribution  of  the  rotation 
frequency  fR. 

The  examples  considered  so  far  have  all  been  numerically  generated. 
Now  examine  the  calculation  of  the  time  domain  average  of  the  vibration  of 
the  Wessex  input  pinion  in  practice  using  sample  and  hold,  linear  and  cubic 
interpolation.  All  examples  are  calculated  using  the  same  input  data 
sampled  and  stored  on  disc  by  a  separate  program.  Samples  were  taken  at  an 
effective  rate  of  20480  Hz  per  channel.  The  vibration  signal  was  lowpass 
filtered  at  4  kHz.  The  results  after  512  averages  are  shown  in  Figure  25, 
26  and  27.  There  is  little  visible  difference.  Previously  a  narrowband 
enhancement  technique  has  been  developed  for  analysis  of  the  vibration  of 
the  Wessex  pinion  [1].  This  technique  can  be  applied  to  averages  presented 
here  by  bandpass  filtering  in  the  range  8  to  36  orders  inclusive  and 
eliminating  the  meshing  component  at  22  orders,  then  enveloping  and 
calculating  the  kurtosis.  The  results  are  shown  in  Table  1.  For  linear 
and  cubic  interpolation,  the  kurtosis  values  obtained  differ  by  only 
0.13  X,  but  the  sample  and  hold  interpolation  differs  from  the  others  by 
more  than  4.9  %.  While  this  difference  is  not  enough  to  affect  the 
decision  on  the  condition  of  the  gear,  it  does  indicate  that  measurable 
changes  in  the  passband  have  occurred.  Views  of  the  spectra  of  the  time 
domain  averages  are  shown  in  Figures  28,  29  and  30,  enlarged  so  as  to 
reveal  the  components  in  the  vicinity  of  the  first  sidelobe.  The  larger 
peaks  in  the  spectra  have  amplitudes  of  approximately  10  g.  For  the  sample 
and  hold  interpolation,  the  components  in  the  sidelobe  have  amplitudes 
about  0.03  g  present  as  a  wideband  noise  across  the  spectrum,  even 
intruding  into  the  passband  and  so  affecting  the  result.  The  best  result 
is  clearly  obtained  by  cubic  interpolation. 
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5.  EXECUTION  TIMES  OF  INTERPOLATION  PROGRAMS 


The  execution  times  of  programs  running  on  a  DEC  LSI  11/73  under  RT-11 
Version  5  with  Fortran  IV  using  sample  and  hold,  linear  and  cubic 
interpolation  to  process  128  averages  of  identical  data  are  tabulated  in 
Table  2.  The  cubic  interpolation  program  requires  more  than  twice  the  time 
needed  for  linear  interpolation.  Table  3  compares  the  execution  times  of 
the  cubic  interpolation  program  to  process  100  averages  of  identical  data 
when  compiled  under  Fortran  IV  and  Fortran  77.  Clearly,  the  greater 
optimization  of  code  performed  by  the  latter  compiler  reduces  the  execution 
time  by  nearly  one  half. 


6.  MANIPULATION  OF  SPECTRUM 


Since  the  time  domain  average  is  exactly  periodic,  the  spectrum  of  the 
average  is  a  pure  line  spectrum.  By  manipulating  components  in  the 
frequency  domain  and  then  transforming  back  to  the  time  domain,  digital 
filtering  can  be  performed  to  improve  the  passband  and  stopband  shapes. 
For  example,  by  setting  the  amplitudes  of  all  the  components  at  frequencies 
in  the  stopband  to  zero,  the  first  sidelobe  can  be  eliminated. 
Furthermore,  if  the  passband  is  divided  by  the  theoretical  frequency 
response  of  the  interpolation  function,  the  passband  can  be  made  flat. 
However,  setting  of  the  stopband  to  zero  will  not  remove  aliasing  errors 
from  the  passband,  and  dividing  the  passband  can  increase  the  amplitude  of 
the  aliasing  errors  there.  Noise  which  is  aliased  into  the  passband  cannot 
be  eliminated  by  manipulation  of  the  passband  shape  as  a  substitute  for  a 
good  interpolation  function. 
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7.  CONCLUSIONS 


Interpolation  techniques  provide  an  alternative  to  the  phase-locked 
frequency  multiplier  for  the  calculation  of  the  time  domain  average  of  a 
vibration  signal  when  a  rotational  reference  signal  taken  directly  from  the 
gear  of  interest  is  not  available.  The  frequency  response  of  the 
interpolation  technique  depends  on  the  order  of  the  technique  used.  Higher 
order  techniques,  such  as  Lagrangian  cubic  interpolation,  have  a  flatter 
passband  and  smaller  sidelobes  in  the  stopband  than  simpler  techniques  such 
as  sample  and  hold  and  linear  interpolation.  Hovever,  increasing  the  order 
of  the  interpolation  technique  also  increases  the  time  required  for  the 
calculation. 

When  the  signal  is  interpolated  at  the  nev  frequency,  the  baseband  of 
the  signal  is  replicated  in  the  frequency  domain  at  multiples  of  the  nev 
frequency.  When  sidelobes  produced  by  the  interpolation  function  are 
replicated  into  the  original  baseband,  aliasing  errors  are  introduced  vhich 
corrupt  the  interpolated  signal.  In  general,  time  domain  averaging 
attenuates  the  aliasing  errors,  but  under  certain  conditions,  it  is 
possible  for  noise  signals  to  generate  aliasing  errors  vhich  vill  be  passed 
by  the  time  domain  averaging  process  vithout  attenuation.  In  practice, 
small  fluctuations  in  the  rotational  speed  of  the  mechanical  system  vill 
assist  the  time  time  domain  averaging  to  attenuate  the  aliasing  errors. 
Nevertheless,  the  use  of  an  interpolation  technique  of  sufficiently  high 
order  is  recommended  to  reduce  the  sidelobes  and  hence  the  aliasing  errors. 
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TABLE  1 


Effect  of  Interpolation  Scheme  on  Enhanced  Signal 
DEC  LSI  11/73  RT-11  Vessex  Input  Pinion  512  averages 


Interpolation  Scheme 

Kurtosis 

Sample  and  Hold 

3.149 

Linear 

3.001 

Cubic 

2.997 

TABLE  2 

Effect  of  Interpolation  Scheme  on  Execution  Time 
DEC  LSI  11/73  RT-11  Fortran  IV  128  averages 


Interpolation  Scheme 

Sample  and  Hold 

Linear 

Cubic 


Execution  Time  (s) 

135 

174 

431 


TABLE  3 

Effect  of  Compiler  on  Execution  Time 
DEC  LSI  11/73  RT-11  Cubic  Intel  )><•!.<  t  .on  100  averages 

Comp l  lei  F.xe, t  ion  Ti iu.  _^s  ) 

Fortran  IV  <  t7 

l  ?4 


Fortran  7 7 


FIG.  2  SAMPLE  AND  HOLD  INTERPOLATION  OF  SAMPLED  SIGNAL 
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FIG.  5  COMPARISON  OF  FREQUENCY  RESPONSES  OF  INTERPOLATION  TECHNIQUES 


PIG  6  SAMPLING  OF  LINEAR  INTERPOLATED  SIGNAL 
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FIG.  7  AMPLITUDE  SPECTRUM  AFTER  SAMPLE  AMD  HOLD  INTERPOLATION  OF 
TEST  SIGNAL  1  AVERAGE 
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FIG.  8  AMPLITUDE  SPECTRUM  AFTER  LINEAR  INTERPOLATION  OF  TEST 
SIGNAL  1  AVERAGE 
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FIG.  9  AMPLITUDE  SPECTRUM  AFTER  CUBIC  INTERPOLATION  OF  TEST  SIGNAL 
1  AVERAGE 


Frequency  (orders) 
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AMPLITUDE  SPECTRUM  AFTER  SAMPLE  AND  HOLD  INTERPOLATION 
OF  TEST  SIGNAL  128  AVERAGES 


FIG.  11  AMPLITUDE  SPECTRUM  AFTER  LINEAR  INTERPOLATION  OF  TEST 
SIGNAL  128  AVERAGES 


FIG.  12  AMPLITUDE  SPECTRUM  AFTER  CUBIC  INTERPOLATION  OF  TEST 
SIGNAL  128  AVERAGES 
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FIG.  13  AMPLITUDE  SPECTRUM  AFTER  SAMPLE  AND  HOLD 
INTERPOLATION  OF  TEST  SIGNAL  1  AVERAGE 


FIG.  14  AMPLITUDE  SPECTRUM  AFTER  SAMPLE  AND  HOLD 

INTERPOLATION  OF  TEST  SIGNAL  123  AVERAGES 


FIG.  15  AMPLITUDE  SPECTRUM  AFTER  LINEAR  INTERPOLATION  OF 
TEST  SIGNAL  1  AVERAGE 


FIG.  1G  AMPLITUDE  SPECTRUM  AFTER  LINEAR  INTERPOLATION  OF 
TEST  SIGNAL  123  AVERAGES 
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FIG.  17  AMPLITUDE  SPECTRUM  AFTER  CUBIC  INTERPOLATION'  OF 
TEST  SIGNAL  1  AVERAGE 
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FIG.  18  AMPLITUDE  SPECTRUM  AFTER  CUBIC  INTERPOLATION  OF 
TEST  SIGNAL  123  AVERAGES 


FIG.  19  AMPLITUDE  SPECTRUM  AFTER  SAMPLE  AND  HOLD  INTERPOLATION  OF 
TEST  SIGNAL  1  AVERAGE 


FIG.  20  AMPLITUDE  SPECTRUM  AFTER  SAMPLE  AND  HOLD  INTERPOLATION 
OF  TEST  SIGNAL  128  AVERAGES 
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FIG.  21  AMPLITUDE  SPECTRUM  AFTER  LINEAR  INTERPOLATION  OF 
TEST  SIGNAL  1  AVERAGE 


FIG.  22  AMPLITUDE  SPECTRUM  AFTER  LINEAR  INTERPOLATION  OF 
TEST  SIGNAL  123  AVERAGES 
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FIG.  23  AMPLITUDE  SPECTRUM  AFTER  CUBIC  INTERPOLATION  OF 
TEST  SIGNAL  1  AVERAGE 


FIG.  24  AMPLITUDE  SPECTRUM  AFTER  CUBIC  INTERPOLATION  OF 
TEST  SIGNAL  123  AVERAGES 
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FIG.  25  TIME  DOMAIN  AVERAGE  OF  WESSEX  INPUT  PINION 

VIBRATION  USING  SAMPLE  AND  HOLD  INTERPOLATION 


FIG.  26  TIME  DOMAIN  AVERAGE  OF  WESSEX  INPUT  PINION 
VIBRATION  USING  LINEAR  INTERPOLATION 
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FIG.  27  TIME  DOMAIN  AVERAGE  OF  WESSEX  INPUT  PINION 
VIBRATION  USING  CUBIC  INTERPOLATION 
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FIG.  28  AMPLITUDE  SPECTRUM  OF  TIME  DOMAIN  AVERAGE  OF  WESSEX 
INPUT  PINION  VIBRATION  USING  SAMPLE  AND  HOLD 
INTERPOLATION 
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FIG.  29  AMPLITUDE  SPECTRUM  OF  TIME  DOMAIN  AVERAGE  OF  WESSEX 
INPUT  PINION  VIBRATION  USING  LINEAR  INTERPOLATION 
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FIG.  30  AMPLITUDE  SPECTRUM  OF  TIME  DOMAIN  AVERAGE  OF 
WESSEX  INPUT  PINION  VIBRATION  USING  CUBIC 
INTERPOLATION 
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